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Abstract 

A theoretical study proposes the atomic configuration of electron-beam irradiated Cgo thin films. 
We examined the electronic structure and electron-transport properties of the Cgo clusters us- 
ing density-functional calculations and found that a rhombohedral Ceo polymer with sp^-bonded 
dumbbell-shaped connections at the molecule junction is a semiconductor with a narrow band gap 
while the polymer changes to exhibit metallic behavior by forming sp^-bonded peanut-shaped con- 
nections. Conductance below the Fermi level increases and the peak of the conductance spectrum 
arising from the t^i states of a Cgo molecule becomes obscure after the connections are rearranged. 
The present rohmbohedral film, including the [2+2] four-membered rings and peanut-shaped con- 
nections, is a candidate to represent the structure of the metallic Ceo polymer at the initial stage 
of electron-beam irradiation. 

PACS numbers: 71.30.+h, 71.15.Mb, 72.80.Rj, 73.61.Wp 



10 I. INTRODUCTION 



11 Intermolecular electron transport has attracted a great deal of attention in numerous fields 

12 of research accompanied by the progress in nanostructure-fabrication technology. Fullerene- 

13 based molecular crystals and films are of particular interest because of their very rich physical 

14 properties such as superconductivity and magnetism. The polymerization of Ceo molecules 

15 in the solid phase using various techniques, e.g., photoirradiation or electron-beam (EB) 

16 irradiation, have resulted in new forms of carbon materials.- The first-principles study by 

17 Okada and Oshiyama reported that the band gap of the rhombhedral phase of the Cqq 

18 polymer linked by [2+2] four-membered rings, where two C atoms shared by the two adjacent 

19 hexagons in a Geo molecule are covalently bonded to the C atoms shared by the two hexagons 

20 in the neighboring Ceo molecule (e.g., see Fig. [1]), is significantly smaller than that of Ceo 

21 bulk.- Onoe et al. found that Ceo films change to polymers accompanying the insulator- 

22 metal transition due to EB irradiation and that a one- dimensional (ID) peanut-shaped 

23 connection is formed, where Ceo molecules are linked by sp^-like connections instead of 

24 sp'^-like connections,- and intensive studies on electron transport in the Ceo polymer have 

25 been carried out up to now.^ First-principles calculations on the electronic structure of 

26 three-dimensional (3D) Ceo polymers have revealed that ID peanut-shaped Ceo polymers 

27 are insulators and electrons are conducted across the ID polymers via [2+2] four-membered 

28 rings between the two peanut-shaped Ceo polymers in the 3D hexagonal structure,-"^ which 

29 is a completely different configuration from the face-centered cubic structure of Ceo bulk. 

30 Nakaya et al. recently revealed that the chemical linkage between Ceo molecules created 

31 by EB irradiation are thermally more stable than the [2+2] four-membered ring because the 

32 polymers are stable under thermal annealing at 220 °C, which is higher than the decomposi- 

33 tion temperature of the [2+2] four-membered rings, and they claimed that a peanut-shaped 

34 fullerene contributes to the metallic characteristics of EB irradiated film.- They also demon- 

35 strated by scanning-tunneling spectroscopy (STS) that the characteristic peak of spectra 

36 of the Ceo molecule disappears at its polymerized center after EB irradiation. Moreover, 

37 their scanning tunneling microscopy (STM) images demonstrated that the Ceo polymers 

38 retain a rhombohedral structure at the initial stage of EB irradiation; however, this result 

39 does not agree with the hexagonal structure derived by Onoe et al. using first-principles 

40 calculations.-i^ Thus, the physics underlying the generation of conductivity in EB irradiated 
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41 Ceo films remains unclear. 

42 In this paper, the atomic configuration of the EB irradiated Ceo film is proposed by using 

43 first-principles calculations. We examine what effect the bond configuration at the molecule 

44 junctions has on the electronic structure and electron-transport properties to interpret the 

45 generation of conductivity in EB irradiated Ceo films. Our findings are that the rhombo- 

46 hedral Ceo polymer consisting of the [2-1-2] four-membered rings and the sp^-like interlayer 

47 connections create a semiconductor with a narrow band gap and the band gap vanishes 

48 when the sp^-like interlayer connection is formed. We found significant differences between 

49 the conduction spectra of the dimers bonded by the sp^-like and sp^-like connections; the 

50 conductance of the sp^-like bonded dimer was higher than that of the sp^-like bonded dimer 

51 below the Fermi level and the peak of the conductance spectrum caused by the t„i orbitals 

52 of the sp^-like bonded dimer was clear while that of the sp^-like bonded dimer was obscure. 

53 All calculations are performed within the framework of density functional theoryS using 

54 the real-space finite-difference approach, which enables us to determine the self-consistent 

55 electronic ground state with a high degree of accuracy by using a timesaving double-grid 

56 technique.-*^ The norm-conserving pseudopotentials^ of Troullier and Martins^^ are used 

57 to describe electron-ion interaction, and exchange correlation effects are treated with local 

58 density approximation.— 



59 II. ELECTRONIC STRUCTURE OF Ceo POLYMERS 

60 Ceo polymers form a triangular lattice with [2-1-2] four-membered rings in each poly- 

61 merized layer and the layers are stacked along a direction perpendicular to the layers in 

62 a rhombohedral symmetry.- When the ^i?C-stacking structure^ is formed, the layers are 

63 bonded by sp^-like connections, in which the hexagons on neighboring Ceo molecules are 

64 brought face to face. When the C-C bonds sheared by a hexagon and pentagon on the 

65 facing hexagons are distorted, the layers are linked by an sp^-like connection consisting 

66 of hexagons and heptagons. After this, we will refer to sp^-like interlayer connections as 

67 dumbbell- shaped and sp^-like interlayer connection as peanut-shaped. The peanut-shaped 

68 connection is much stronger than the [2+2] four-membered ring because the energy to form 

69 it is smaller than that for the [2-|-2] four-membered ring by 0.02 eV/atom.^^ Therefore, the 

70 polymer can be reinforced by the peanut-shaped connection so that it can survive under 
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71 thermal annealing. We first calculate the electronic structure of the Ceo polymers. 

72 Figure [T] shows two geometries of the models studied here. The Cgo polymers in model (a) 

73 are bonded by the [2+2] four-membered rings in the layer and linked by the dumbbell-shaped 

74 connections between the layers. In model (b), one of the dumbbell- shaped connections 

75 between the layers is deformed into the peanut-shaped connection. Since the metallic phase 

76 of the Ceo polymer is surrounded by the other phase of the Ceos at the initial stage of 

77 the EB irradiation according to the STM image/ its lattice parameters are not expected 

78 to be fully relaxed. We assume lattice parameters a=17.37 bohr, c=46.30 bohr and a = 

79 7r/3, which correspond to those for the dumbbell- shaped interlayer connection obtained 

80 by the x-ray diffraction pattern analysis,— and employ these lattice parameters for both 

81 model (a) and model (b) to compare the contribution of the dumbbell-shaped and peanut- 

82 shaped connections to the electronic structures. Note that the stacking structure of the 

83 polymer presented in this study is not deformed from the rhombohedral bulk phase and 

84 corresponds to the experimentally observed STM image whereas the hexagonal polymers in 

85 Refs. al and require drastic structural deformation. Integration in the Brillouin zone is 

86 carried out using 24-point sampling and structural optimization for the atomic geometris 

87 is implemented for both models until the remaining forces for each atom are less than 36 

88 meV/bohr. The calculated electronic band structures are depicted in Fig. O The energy 

89 band gap is reduced by forming peanut-shaped connections; the peanut-shaped polymer 

90 exhibits metallic behavior while model (a) is a semiconductor with a fundamental band gap 

91 of ~ 0.6 eV. The formation energy of model (b) is larger than that of model (a) by 0.25 

92 eV/atom although the formation energy of an isolated peanut-shaped dimer is smaller than 

93 that of a dumbbell-shaped dimer by 0.04 eV/atom. Thus, the peanut-shaped connections are 

94 not widely generated but Ceo molecules form small clusters with peanut-shaped connections 

95 in the film after EB irradiation. 



96 III. TRANSPORT PROPERTIES OF Cgo DIMERS 

97 The STS spectra, which have been employed to interpret the transport properties of 

98 molecules, show peaks which are distinctly associated with the electronic structure of 

99 molecules.— However, to replicate the STS spectrum of the Cqo cluster in the Cgo film, the 

100 system including the large number of atoms are required because a Ceo molecule contains 



101 60 atoms, which is not easy task by the present computational resources. Moreover, the 

102 number of peaks measured for the considered bias-voltage range have been smaller than 

103 the number of electronic states in the molecules within the corresponding energy window in 

104 some cases.— This fact implies that some states of molecules do not contribute to electron 

105 transport and that a simple interpretation of transport properties only in terms of DOS 

106 is insufficient. In addition, it is well known that the junction between Ceo molecules are 

107 a bottleneck for electron transportp^"— and the energy gap between the highest-occupied 

108 molecular orbital (HOMO) and the lowest-unoccupied molecular orbital (LUMO) decreases 

109 after polymerization.-i^ii^i^ One of the present authors (S.T.) examined variations in the 
no local density of states (LDOS) along the ID Cqo polymer axis and demonstrated that the 

111 energy gap of LDOS at the molecule junction is larger than that in the Cqo molecules.— 

112 The sp^-like connections of the peanut-shaped polymer might enhance electron transport 

113 because vr electrons are major carriers in metallic carbon nanotubes. Thus, it is of con- 

114 siderable interest to examine how the difference in the bond network between fuUerenes 

115 contributes to electron transport through fuUerenes. 

116 We explore the contribution of the peanut-shaped connection to electron-transport prop- 

117 erties. Figure [3] shows the computational model, where a Cgo dimer is sandwiched between 

118 electrodes. Since the lattice-constant mismatch between the Ceo polymers and the elec- 

119 trodes gives rise to further complex discussions on the transport properties, we employ the 

120 simplified models to focus on the difference in the contributions to the transport properties 

121 between the dumbbell-shaped and peanut-shaped connections. To determine the optimized 

122 atomic coordinates and Kohn-Sham effective potential, we use a conventional supercell un- 

123 der a periodic boundary condition in all directions with a real-space grid spacing of ~ 0.32 

124 bohr; the dimensions of the supercell are = 37.87 bohr, Ly = 37.49 bohr, and Lz = 71.15 

125 bohr, where and Ly are the lateral lengths of the supercell in the x- and y-directions 

126 parallel to the electrode surfaces, respectively, and is the length in the z direction. 

127 Structural optimizations are implemented in advance for the isolated peanut-shaped dimer 

128 and it is then placed between the electrodes, where the three topmost surface atomic lay- 

129 ers are atomistic Al(lll) and the rest are aluminum jellium. The dimer is aligned at the 

130 hexagonal-close-packed-hollow site on the (111) surface facing a hexagon, which is the most 

131 stable configuration for a Ceo molecule on a face centered cubic (111) surface according to 

132 first-principles calculations.— The distance between the surface atomic layers on the left 
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133 and right electrodes is set at 35.81 bohr so that the distance between the edge atoms of the 

134 dimer and the surface atomic layer of the electrodes corresponds to that reported by first- 

135 principles calculations. The dimer is relaxed between the electrodes. The atomic geometries 

136 of 36 carbon atoms for the dumbbell-shaped dimer at the molecule junction are modified 

137 from the above mentioned structure and structural optimization is also implemented. 

138 We take a grid spacing of ~ 0.47 a.u for the electron-transport calculations. We ensured 

139 that the decreased grid spacing and the enlarged supercell would not significantly affect our 

140 results. The scattering wave functions of the electrons propagating from the electrodes are 

141 determined by using the method of overbridging boundary matching.—"— The retarded self- 

142 energy matrices for aluminum jellium are employed to include the rest of the semi-infinite 

143 electrodes. Since the DOS of aluminum is similar to that of free electrons, unfavorable 

144 effects from the DOS of the electrodes on the conductance spectra can be eliminated. We 

145 first calculate the Kohn-Sham effective potential using the supercell employed in structural 

146 optimization and then compute the scattering wave functions obtained non-self-consistently. 

147 It has been reported that this procedure is just as accurate in the linear response regime but 

148 significantly more efficient than performing computations self-consistently on a scattering- 

149 wave basis.— The conductance of the dimers is described by the Landauer-Biittiker formula, 

150 G = Tr(T^T)G'o;™ where T is a transmission-coefficient matrix. 

151 Figure m plots the conductance spectra of the dumbbell-shaped and peanut-shaped dimers 

152 as a function of the energy of incident electrons. The magnitude of conductance at the Fermi 

153 level and the conduction spectra as a function of the energy of the incident electrons for 

154 the dumbbell-shaped dimer are in agreement with those reported in another theoretical 

155 study.— Since the dimers are not connected by the [2+2] four-membered rings in a direction 

156 parallel to the electrode and the ID Ceo polymers act as an insulator, the conductance 

157 of the dimers is low around the Fermi level. This implies that the [2+2] four-membered 

158 rings play an important role in the metallic characteristics of the Cgo polymers and the 

159 very bright Cgo molecule in the experimental STM image- is presumably linked by the 

160 [2+2] four-membered rings and the peanut-shaped connections as explained in the preceding 

161 section. The difference in energy between the peaks of the spectrum around the Fermi level 

162 corresponds to the energy gap between the HOMO and the LUMO of the isolated dimers. 

163 The peak at Ep + 0.25 eV is attributed to the triply degenerate LUMO (tui) of a Ceo 

164 molecule, where Ep is the Fermi level. The conductance of the peanut-shaped dimer is 
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165 higher than that of the dumbbell-shaped dimer below the Fermi level and the peak of the 

166 conductance spectrum induced by the t^i orbitals of the dumbbell-shaped dimer is clearer 

167 than that of the peanut-shaped dimer. 

168 To examine the contribution of the peanut-shaped connection to the electron-transport 

169 properties in more detail, Fig. [5] shows the LDOS of the dimers, which have been plotted by 

170 integrating them along the plane parallel to the electrode surface, p{z, E) = J \ip{r, E)\'^dr\\, 

171 where r = {x, y,z), ip is the wave function and E is the energy of the states. The LDOS at 

172 Ep — 0.2 eV around the molecule junction of the peanut-shaped dimer is larger than that 

173 of the dumbbell-shaped dimer. The energetically discrete t„i orbitals of the Ceo molecules 

174 contributing the electron transport change to the broadened states because of the formation 

175 of the sp^-like connection at the molecule junction. The states in the electrode can easily 

176 penetrate into the peanut-shaped dimer and the conductance spectrum in Fig. H] becomes ob- 

177 scure. Although our computational models do not directly correspond to the tip and surface 

178 system of STS, the absence of clear peaks of the tui orbitals at the center of the polymerized 

179 cluster in the STS spectrum is related to the smooth behavior of the conductance spectrum 

180 induced by the sp^-like peanut-shaped connections. 



181 IV. SUMMARY 

182 The atomic configuration of the EB irradiated Ceo polymers has been proposed to ex- 

183 plain the generation of the conductivity in the Ceo fil^i through first-principles calculations. 

184 Although the polymer composed of the [2-1-2] four-membered rings and dumbbell-shaped 

185 interlayer connections is a semiconductor with a narrow band gap, the polymer changes to 

186 exhibit metallic characteristics by forming a peanut-shaped interlayer connection. The cal- 

187 culations for electron transport for the dumbbell-shaped and peanut-shaped dimers revealed 

188 that the higher conductance of the peanut-shaped dimer below the Fermi level than that of 

189 the dumbbell-shaped dimer and the low peak of the conductance spectrum induced by the 

190 tui orbitals of the peanut-shaped dimer are relevant to the energetically broadened LDOS at 

191 the molecule junction caused by the formation of the sp^-bonded peanut-shaped connection 

192 and the metallic property of the EB-irradiated Ceo polymer. Rhombohedral film consisting 

193 of [2-1-2] four-membered rings, dumbbell-shaped and peanut-shaped connections is a possible 

194 structure for the metallic Ceo polymer observed by STM at the initial stage of EB irradia- 
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195 tion, since the present structure could be formed without significant structural deformation 

196 from face-centered cubic Ceo bulk and agrees with the rhombohedral structure in the STM 

197 image. Although our computational models are limited to the simplified Ceo clusters because 

198 of the limitation of the computational resources, the large-scale first-principles calculations 

199 for the Ceo polymers could help to validate the correlation between the STS spectra and the 

200 present atomic configuration of the polymers and will be carried out in a future. 
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FIG. 1. (color online) Geometric structure of 3D Cqo polymer phases, (a) Polymer bonded by only 
[2+2] four-membered rings and dumbbell-shaped connections and (b) polymer in which one of three 
dumbbell-shaped connections deforms to the peanut-shaped connection. The shaded segments in 
(b) represent the hexagons (blue) and heptagons (red) of the peanut-shaped connection. B and G 
in the top views correspond to the positions of Ggo molecules in the upper and lower layers. 
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FIG. 2. Band structures of 3D Cqo polymers of models (a) and (b) in Fig. [TJ Zero energy was 
chosen as the Fermi level. 
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FIG. 3. (color online) Computational model where Cqo dimer is suspended between Al(lll) elec- 
trodes, (a) Ddumbbell-shaped dimer and (b) peanut-shaped dimer. 
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FIG. 4. Conductance spectra as function of energy of incident electrons. The dashed curve repre- 
sents the conductance of the dumbbell-shaped dimer and the solid curve that of the peanut-shaped 
dimer. Zero energy was chosen as the Fermi level. 
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FIG. 5. (color online) Distributions of LDOS integrated on plane parallel to dimer as functions of 
relative energy from Fermi level. Zero energy is chosen as the Fermi level. Each contour represents 
twice or half the density of the adjacent contour lines, and the lowest contour is 2.56 xlO~^ 
e/eV/bohr. The atomic configurations below the graph are visual guides. 
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